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PREFACE 


As  a  pare  of  RAND's  continuing  interest  in  weather  modification 
in  its  broadest  aspects,  the  construction  and  study  of  numerical  models 
of  atmospheric  orocasees  on  several  scales  has  been  emphasized.  An  im¬ 
portant  part  of  this  activity  has  been  the  develop.nent  of  a  model  treat¬ 
ing  the  hydrodynamics  and  thermodynamics  of  an  individual  cumulus  cell. 

A  number  of  investigators  have  addressed  themselves  to  this  prob¬ 
lem,  and  most  have  been  constrained  by  computer  capacity  to  assume  some 
form  of  symmetry  in  order  to  reduce  the  problem  from  three  dimensions 
to  two.  It  has  generally  been  considered  that  formula'  an  in  Cartesian 
coordinates  is  satisfactory,  although  the  more  complex  cylindrical  for¬ 
mulation  more  nearly  describes  what  occurs  in  nature.  Except  for  a 
brief  discussion  by  Ogura  in  1963,  however,  the  two  formulations  have 
not  been  directly  compared.  It  is  the  purpose  of  this  study  to  make 
such  a  comparison  in  greater  depth  and  so  to  determine  the  conditions 
under  which  each  of  the  formulations  is  best  used. 

Other  RAND  publications  related  to  the  present  study  are  RM-5316-NRL, 
Numerical  Simulation  of  Cumulus  Convection;  RM-5564-NRL,  A  Method  of 
Objective  Contour  Construction;  and  RM-5582-ESSA ,  An  Annotated  Bibli- 


ABSTRACT 


Two  versions  of  a  numerical  model  for  cumulus  convection  are  com¬ 
pared.  One  is  symmetrical  about  a  vertical  plane*  the  other  about  a 
vertical  axis.  It  is  found  that  the  axisymmetr ic  model  glows  more  vig- 
orously  than  the  other,  but  more  realistically  represents  the  relations 
between  updraft  an^  downdraft,  the  shape,  and  other  characteristics . 

The  model,  at  least  with  the  eddy  exchange  coefficients  used,  appears 
to  be  deficient  in  turbulent  entrainment.  The  previous  findings  of 
Ogura  are  generally  confirmed  and  extended. 


I.  INTRODUCTION 


Since  the  initial  nume'ical  experiments  in  atmospheric  thermal 
convection  by  Malkus  and  Witt  (1959)  and  Ogura  (1962,  1963)  a  number 
of  otner  models  of  the  same  general  type  have  been  proposed,  each  em¬ 
phasizing  one  or  another  aspect  of  the  problem.  Among  theae  are  the 
models  of  Ogura  and  Charney  (1962),  Lilly  (1962),  Chou  (1962),  Chou, 

Li,  Chang,  and  Kung  (1964),  Li,  Chao,  and  Hwu  (  1964),  Asai  (1964a, b', 
Murray  and  Anderson  (1965),  Nickerson  (1965),  Takeda  (1965),  Orville 
(1965  1968),  Amirov  (1966a, b),  Inman  (1966),  Lebedev  (1966),  and  'Arns 

son,  Greenfield,  and  Newburg  (1968).  These  models  differ  in  many  de¬ 
tails,  but  nearly  all  of  them  have  in  common  the  property  of  two-di¬ 
mensionality,  the  full  extension  to  th’-ee  dimensions  being  beyond  the 
capacity  of  computers  available  during  this  period. 

There  have  been  two  principal  v.ay.  in  which  the  third  dimension 
was  eliminated.  In  some  cases  ('artesian  coordinates  were  used  with 
the  assumption  that,  no  property  varies  in  the  y-dlrection,  whereas  in 
others  cylindrical  coordinates  were  used  ~ii.ii  the  assumption  that  no 
property  varies  in  the  azimuthal  direction.  Certain  problems  ( e .  g ,  , 
the  squall  lino)  virtually  demand  the  ('artesian  approach  ,  hut  for  model 
fng  individual  clouds  the  cylindrical  formulation  is  more  ippropr i ate . 
However,  the  additional  complexity  of  the  vlindrical  formulation  has 
led  some  investigators  to  use  the  other,  even  when  it  Is  apparently 
less  suited  to  the  problem  at  hand.  So  far  as  has  been  ascertained, 
only  Ogura  (19t>J)  has  studied  the  same  problem  wl  th  both  georx'tr  les , 
and  it  is  «'ii  the  basis  of  rather  limited  comparative  data  published  bv 
him  that  it  has  been  assu.ied  that  the  two  schemes  yield  essentially 
similar  results.  It  is  the  purpose  oi  the  present  study  to  farther 
test  this  issumptlon,  and  to  do  so  it  is  desirable  to  make  a  more  de¬ 
tailed  comparison  of  axial  and  rectilinear  symmetry  in  numerical  cloud 
mode  ling. 

Most  of  the  models  mentioned  above  are  based  on  the  Bouasinesq  ap¬ 
proximation  or  something  very  close  u  it,  which  permits  the  definition 
of  a  stream  function  and  the  derivation  of  a  vorticitv  equation.  A 
number  of  1  ini te -di i t e rence  schemes  have  been  devised  to  solve  this 
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type  of  equation.  The  properties  of  thee  schemes  have  been  extensively 
discussed  in  the  literature  (see,  e.g.,  Molencamp,  1968)  and  will  not 
be  further  discussed  here.  Almost  all  of  the  prr-ticai  schemes,  how¬ 
ever,  solve  the  vorticity  equation  locally  in  the  Eulerian  sense,  fol¬ 
lowing  which  the  stream  function  is  found  by  inversion  of  an  elliptic 
operator,  usually  by  an  iterative  technique. 

Several  o*  the  convection  models  treat  a  dry  atmosphere.  This  ap¬ 
proach  can  yield  some  interesting  results,  but  it  is  much  more  useful 
to  consider  a  moist  atmosphere  in  which  water  can  exist  in  at  least  the 
liquid  and  vapor  phases.  Murray  and  Hollinden  (1966)  have  proposed  a 
model  that  treats  water  in  all  three  phases.  It  has  been  found  that 
in  atmospheric  thermal  convention,  release  of  latent  heat  has  more  ef¬ 
fect  on  buoyancy  than  any  other  single  process.  Some  models  assume 
pseudoad labatic  processes,  i.e.,  all  condensed  water  is  immediately 
dropped  out.  Hence  a  cloud  parcel  ascends  moist-adiabatical ly ,  but 
descends  dry-adlabaticr.l ly  without  evaporation.  Other  models,  includ¬ 
ing  the  one  to  be  described  herein,  assume  reversible  processes;  i.e., 
all  condensed  water  is  carried  with  the  air  and  is  available  later  for 
evaporation.  Still  better  is  the  assumption  that  only  part  of  the  con¬ 
densate  falls  out  as  rain,  but  determination  of  what  part  this  is  to 
be  is  still  a  vexing  problem.  Precise  treatment  of  rain  would  necessi¬ 
tate  detailed  modeling  of  mfcrophysical  processes,  but  to  date  no  pro¬ 
gram  has  been  developed  that.  can  handle  both  this  and  cloud-scale  dy¬ 
namics  simultaneously :  in  all  probability  this  must  await  the  next  gen¬ 
eration  of  large,  psral le 1 -processing  computers.  Several  existing 
models,  however ,  have  parameterized  precipi tation  mechanisms.  The  most 
successful  parameterization  of  this  type  so  far  is  the  one  developed  1. 
Kessler  (1967).  Hie  model  described  herein  is  based  on  reversible  pro¬ 
cesses  without  fallout  of  rain. 

The  method  developed  by  Ogura  (  1963 '  to  handle  moist  air  (and  later 
»u.  ted  by  several  other  investigators  1  is  is  follows:  From  the  basic 
equations  of  hydrodynamics  and  thermodynamics,  one  derives  a  prediction 
equation  for  the  meridional  component  of  vorticity  and  conservation 
equations  for  total  water  substance  and  specific  entropy  of  moist  air. 
The  finite -di f ference  analogs  of  these  equations  are  solved  at  each 
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grici  point,  and  the  temperature  is  then  determined  from  the  specific 
entropy  by  one  of  two  irvethods,  depenuirtg  on  whether  or  not  the  air  is 
saturated.  Total  >.uter  substance  i-  then  partitioned  between  liquid 
and  vapor  on  the  assumption  that  any  excess  of  mixing  ratio  over  the 
saturation  value  represents  liquid  water. 

In  the  present  investigation  a  somewhat  different  method  was  em¬ 
ployed.  The  vorticity  equation  and  its  solution  are  essentially  the 
same  as  in  Ogura's  model,  but  the  thermodynamic  equations  are  solved 
in  a  l.agrangian  manner.  The  rationale  for  this  is  that  the  physical 
processes  relating  to  condensation  ana  evaporation  are  usually  described 
in  terms  of  changes  occurring  j.n  individual  parcels.  The  various  pro¬ 
cesses  are  interrelated  in  a  highly  implicit  way,  but  if  attention  is 
confined  to  an  individual  parcel,  they  can  be  treated  serially  while 
maintaining  the  erro:  within  acceptable  bounds .  The  l.agrangian  approach 
is  conceptually  more  straightforward  than  the  tulerian  in  this  •  ase  , 
and  in  practice  it  has  been  found  to  be  more  stable  computationally. 

One  drawback  is  that  the  smoothing  resulting  from  the  necessary  ‘  iter- 
polation  acts  'ike  »n  impli.it  eddy  d i f i us . on .  This  effect,  however, 
appears  to  ue  snail  enough  not  to  be  troublesome .  Quantities  such  as 
total  water  substance  are  not  exactly  conserved,  but  the  ios*.  is  quite 
sma 1 1  . 
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II.  THE  MODEL 


The  Boussinesq  equation  of  motion  for  moist  air  (valid  if  convec¬ 
tion  is  confined  to  a  layer  a  few  kilometers  in  depth)  is 


( 1 ) 


The  symbols  are  defined  in  the  Appendix.  The  virtual  temperature 
is  defined  by 


* 
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The  Boussinesq  equation  of  continuity  is 


V-V  =  0 

If  one  takes  the  curl  of  (1)  and  uses  (3),  there  results 


(2) 


(3) 
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B  =  g 
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(4) 


(5) 


(6) 


Two  different  .wo-dimensiona 1  models  derived  from  (4)  will  be  dis- 

The  first:  uses  a  Cartesian  coordinate  system  (x,y  ,z )  and  proceeds 


cussed . 
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on  the  assumption  that  the  component  of  velocity  and  the  variation  of 
all  properties  in  the  y-direction  are  zero',  this  is  the  rectilinear 

-u 

model.  Hie  second  uses  a  cylindrical  coordinate  system  (>:,2,z)  and 
proceeds  on  the  assumption  that  the  tangential  component  of  velocity 
is  zero  and  all  properties  are  indeoendent  of  S;  this  is  the  axisyffi- 
metric  model.  Subscripts  R  and  A  will  be  used  to  denote  variables  per¬ 
taining  to  the  rectilinear  and  axisymmetric  models,  respectively. 

Because  of  the  incompressibility  implied  by  (3)  it  is  possible  to 
express  the  wind  components  in  terms  of  a  stream  function,  as  follows: 


(V) 


The  horizontal  component  of  vorcicity  is 


k 

The  coordinate,  x.  is  used  here  rather  than  the  more  conventional 
r  to  facilitate  comparison  with  the  Cartesian  system. 
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The  transverse  or  tangential  component  of  (4)  is  then 


(9) 

If  at  a  given  time  ijf  (which  determines  u,  w,  and  Tj)  and  T,  r  , 
and  rf  (which  collectively  determine  B)  are  known  at  each  p>:id  point, 

(9)  can  be  solved  to  determine  1]  at  each  grid  point  at  a  later  time,  ! 

following  which  (Sa)  can  be  solved  for  V  at  the  later  time.  The  latter  j 

can  be  done  as  a  Dirichlet  problem  with  V  =  0  on  all  boundaries,  a  con¬ 
dition  equivalent  to  the  requirement.  that  no  air  pass  through  any  of 
the  boundaries,  (The  boundaries  are  taken  to  be  two  horizontals,  one  j 

representing  the  ground  and  one  an  undisturbed  level,  and  two  verticals,  ] 

one  representing  an  axis  of  symmetry  and  one  a  lateral  wall  enclosing  i 

Che  region  of  convective  circulation,)  ■ 

The  Lagrangian  part  of  the  scheme  encompasses  the  determination  j 

of  T,  rv,  and  rf  at  each  grid  point  at  the  later  time.  This  is  done 

in  several  atepe,  as  follows: 

1,  Using  the  values  of  u  and  w  at  a  grid  point  3nd  its  neighbors  i 

at  the  beginning  and  end  of  the  time  step,  a  single-time^step  backward  ' 

trajectory  is  computed  to  determine  the  point:  of  origin  cf  the  aii  par-  « 

cel  that  is  found  at  the  grid  point  at  the  end  of  the  time  step, 

2,  Values  of  temperature,  humidity,  etc,  at  the  point  of  origin 
at  the  beginning  of  the  time  step  are  found  by  interpolation  among  the 
four  surrounding  grid  points. 

J.  If  the  air  parcel  at  the  point  of  origin  is  subsaturated  in 
the  presence  of  liquid  water  or  is  supersaturated,  an  adjustment  is 
made  to  T,  and  to  rectify  this  condition.  Since  the  adjustment 
is  swell,  it  is  assumed  to  occur  Instantaneously  at  constant  pressure. 

"aturation  mixing  ratio  is  given  b 
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rs  p  -  e 


and  the  saturation  vapor  pressure  by  Tetens'  formula: 


ra(T  -  273.16: 


es  =  6 « 1078  exPL  (TTit 


with  a  =  17.2693882  and  b  =  35.86.  Differentiation  and  combination  of 
(10)  and  (11)  yield 


r  (e  +  r  )a(273.16  -  b) 

dr  =  — - - - ^ - -  dT 

S  € (T  -  b>2 


But  the  change  of  temperature  due  to  condensation  or  evaporation  is 


dT  =■  - - - - - dr  *  -  - — - r - dr 

c+rc  +  r  „c  l  c+rc  +r.c  v 
p  v  pv  tv  p  v  pv  l  w 


The  object  of  the  adjustment  is  to  make  rg  and  r  equal;  hence 


r  +  dr  *  r  +  dr 
s  s  v  v 


Combination  of  (12),  (13),  and  (14)  yie’ds 


r  -  r 

.  _  .  _ v  s  _ 

t  ~  v  '  r  (e  +  r  )a(273.16  -  b)L 

1  +  ~  - 

«(T  -  b)  <cp  +  Vpv  *r,cw) 


The  static  adjustment  consists  of  applying  (15/  to  r  and  r  and  (13) 

I  V 

to  T.  If  r^  <  rg,  the  adjustment  can  be  no  greater  than  r^,  (.15;  not¬ 
withstanding. 

4.  If  the  parcel  is  now  exactly  saturated,  a  dynamic  adjustment 

is  made  to  T,  r  ,  and  r..  Under  the  assumption  of  conservation  of 
v  l 
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equivalent  potential  temperature,  it  can  be  shown  that  the  rate  of 
change  of  saturation  mixing  ratio  is 


c  T  -  Lr 
1  -  -E - 1 

eL 


gw 


L  + 


c  RT 

±. 


Lr  (e  +  r  ) 
s  s ' 


(16) 


If  r  *  r  ,  then  dr  /dt  *  -dr  /dt  =  dr„/dt,  and  so  (16)  gives  the  rate 
vs  s  v  l 

of  condensation.  Thus  one  can  determine  the  change  of  T,  r^,  and 
due  to  condensation  or  evaporation  as  the  parcel  moves  from  its  point 
of  origin  to  the  grid  point. 

5.  The  temperature  of  the  parcel  also  changes  adiabatica lly  at 
the  rate  -gw/c  ;  this  must  be  taken  into  account  whether  or  not  step 

p 

(4)  was  applicable.  At  this  stage  the  values  of  T,  r  ,  and  r  result- 

V  X 

ing  from  the  application  of  the  various  adjustments  to  the  values  de¬ 
termined  for  the  point  of  origin  at  the  beginning  of  the  time  step  are 
preliminary  values  for  the  grid  point  at  the  end  of  the  time  step. 

2 

6.  Fickian  eddy  diffusion  is  assumed,  so  terms  of  the  form  v  V  T  , 

2  2  1 
r^,  and  v^V  must  be  included  in  the  rates  of  changes  of  tempera¬ 
ture,  mixir.3  ratio  of  vapor,  and  mixing  ratio  of  liquid.  These  terms 
are  evaluated  in  the  Eulerlan  manner  at  each  grid  point  at  the  beginning 
of  the  time  step,  and  are  now  used  co  adjure  the  grid-point  values  of 

r  ,  and  r.  at  the  end  of  the  time  step. 

V  K 

7.  Because  of  nonlinearities  in  most  of  tKe  processes,  there  may 
at  this  stage  be  points  at  which  the  air  is  subsaturated  In  the  pres¬ 
ence  of  liquid  water  or  is  supei saturated.  If  so  the  adjustment  of 
step  (3)  is  again  applied,  this  time  at  f  •«  grid  point. 
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III.  INITIAL  CONDITIONS 

It  is  assumed  that  initially  the  atmosphere  is  at  rest  and  hori¬ 
zontally  homogeneous.  Ordinarily  the  variation  of  T  and  r  in  the  ver¬ 
tical  is  in  accordance  with  seme  actual  atmospheric  sounding,  and 
is  taken  to  be  zero  everywhere.  Unless  some  impulse  is  applied,  the 
right-hand  side  of  (9)  is  identically  zero,  and  nothing  can  happen. 

The  usual  impulse  used  for  this  type  of  model  is  a  bubble  of  increased 
temperature.  This  makes  3  no  longer  independent  of  x,  so  11  can  grow. 

The  rising  motion  so  induced  in  the  center  of  the  bubble  can  lead  to 
condensation  and  conversion  of  latent  heat  to  sensible  heat,  making 
the  piocess  self-sustaining.  Something  like  this  process  is  thought 
to  occur  in  nature  when  a  cumulus  cloud  is  generated,  but  in  one  re¬ 
spect  this  process  is  unrealistic.  The  arbitrary  initial  increase  in 
temperature  within  the  bubble  is  accompanied  by  a  decrease  in  relative 
humidity,  delaying  the  onset  of  condensation  and  ensuring  that  the  cloud 
develops  on  the  driest  part  of  the  region. 

It  is  perhaps  more  realistic  to  introduce  a  humidity  impulse  rather 
than  a  temperature  Impulse.  There  Is  then  a  horizontal  gradient  of 
virtual  temperature  (though  not  of  temperature)  to  make  the  right-hand 
side  of  (9)  nonzero.  If  the  impulse  is  such  as  to  produce  saturation 
within  the  bubble,  condensation  can  start  immediately.  Humidity  im¬ 
pulses  have  been  studied  both  theoretically  and  observatlonally  by 
Vul'fson  (1963),  who  found  them  to  be  of  great  importance,  especially 
in  convection  over  a  water  surface.  This  is  the  type  of  impulse  used 
in  the  present  study. 
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IV.  DISCUSSION 


Ogura  (1963)  presented  results  of  several  runs  with  a  model  having 
axial  symmetry.  He  used  a  30-by-3Q  grid,  and  for  lateral  boundary  con¬ 
ditions  he  required  that  both  u  3  0  and  dw/dx  =  0  at  the  central  axis 
and  at  the  outer  wall.  One  run  having  a  basic  atmosphere  with  a  con¬ 
ditionally  unstable!  lapse  rate  and  saturation  throughout  was  repeated 
with  rectilinear  geometry,  but  with  a  20-by-20  grid  and  cyclic  lateral 
bound,  ry  conditions. 

The  main  points  noted  by  Ogura  in  comparing  the  two  runs  were : 

a.  The  strength  of  the  downdraft  outside  the  cloud  is  stronger 
in  the  rectilinear  model  than  In  the  axisynraetric .  At  8  minutes  the 
ratios  of  maximum  downdraft  to  maximum  updraft  were  .25  and  .06,  re¬ 
spectively. 

b.  The  heating  of  the  subsiding  air  produced  a  temperature  grad¬ 
ient  such  that  a  cell  with  reversed  circulation  appeared  below  the 
cloud  In  the  rectilinear  model,  but  not  In  the  axlsymmetr ic . 

c.  The  greater  excess  temperature  produced  in  the  downdraft  of 
Che  rectilinear  model  resulted  In  appearance  of  a  region  of  positive 
vorti.».ity.  This  did  not  occur  with  the  axisymnetric  model. 

It  Is  possible  that  Ogura's  use  of  different  boundary  conditions 
and  grid  dimensions  contributed  to  the  different  results  noted  in  the 
models  with  the  two  geometries.  In  order  to  avoid  this  possibility  in 
the  present  investigation,  runs  were  made  with  the  two  geometries  using 
identical  Initial  data,  impulses,  boundary  conditions,  and  grid  dimen¬ 
sions.  The  Initial  sounding  used  was  one  from  San  Juan,  Puerto  Rico, 
for  a  day  on  which  well-developed  cumuli  were  observed  in  the  tegion. 

A  humidity  Impulse  was  chosen  so  that  initially  there  was  slight  super- 
saturation  near  the  central  axis  and  the  convective  condensation  level; 
hence  condensation  commenced  iraned .acely,  and  the  model  cloud  grew  vig¬ 
orously. 

The  moat  obvious  difference  in  the  response  of  the  axlsymmetr ic 
and  rectilinear  models  is  that  the  former  developed  much  faster  and 
more  strongly.  This  is  apparent  from  Fig.  1,  which  shows  the  magni¬ 
tude  of  the  maximum  updraft  as  a  function  of  time.  By  8  minutes  the 
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axisymnetr ic  updraft  was  twice  as  strong  as  the  rectilinear,  and  it 
continued  to  increase  exponentially  until  computation  was  stopped  at 
26  minutes,  when  the  updraft  was  nearly  12  m  sec  \  By  contrast,  the 
updraft  nf  the  rectilinear  model  did  not  quite  reach  2  m  sec  ^  in  the 
same  time,  and  its  rate  of  growth  showed  signs  of  leveling  off  at  29 
minutes. 

This  behavior  could  be  predicted  qualitatively  from  the  nature 
of  the  two  geometries.  Under  the  assumption  that  the  horizontal  area 
within  radial  distance  X  is  nX“  for  the  axisymmetric  model  and  2XY  ' o. 
the  rectilinear  model,  where  Y  is  some  distance  along  the  v-axis,  one 
can  determine  that  if  there  is  updraft  from  the  axis  out  to  and  there¬ 
after  downdraft  out  to  X2 ,  the  proportion  of  the  total  horizontal  area 
occupied  by  updraft  13  X^/X^  for  the  axisymmetric  model  and  X,/X,  for 
the  rectilinear.  In  the  present  experiment,  X,  ~  3800  m  for  both  cases, 
whereas  X^  is  typically  800  m  for  the  axisymmetric  case  and  600  m  for 
the  rectilinear.  Thus  the  ppdraft  covers  2  percent  of  the  total  area 
in  the  axisymnetr ic  case  and  10  percent  in  the  rectilinear.  Obviously, 
if  the  two  updrafts  are  to  transport  similar  amounts  of  air  through  a 
given  level,  that  of  the  axisymmetric  model  must  be  considerably  stronger 
than  the  other.  It  should  bt.  remarked,  however,  that  In  both  models 
moat  of  the  downdraft  is  concentrated  near  the  evaporating  edge  of  the 
cloud,  and  the  remainder  of  the  downdraft  area  has  a  very  slight  ver¬ 
tical  motion,  and  even  some  pockets  of  weak  updraft. 

In  his  initial  exposition  of  the  slice  method,  Bjerknes  (1938'' 
showed  that  the  ratio  of  updraft  to  downdraft  area  must  be  small,  and 
in  an  elaboration  of  that  work,  Fetter ssen  (193Q)  suggested  that  10 
percent  might  be  a  reasonable  figure,  although  he  did  not  express  high 
confidence  in  that  estimate.  However,  it  agrees  with  the  present  re¬ 
sults  for  the  rectilinear  model.  Recently  Krishnamurti  (1963)  deter¬ 
mined  by  a  study  of  total  moisture  convergence  over  an  area  of  synoptic 
scale  in  the  Caribbean  and  by  airborne  radar  observations  that  1  per¬ 
cent  of  the  total  area  was  covered  by  active  convective  cloud  elements. 

Hi  is  agrees  with  results  of  the  axisymmetric  model  in  the  present  ex¬ 
periment  and  offers  a  degree  of  confirmation  of  the  superiorly  of  thi> 
model  in  this  respect. 
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One  difficulty  that  occurs  throughout  this  problem  is  the  mathe¬ 
matical  singularity  on  the  central  axis  in  the  axisymmetric  model. 

This  has  caused  the  greatest  trouble  in  the  evaluation  of  w  from  (7) 
where  x  *  0,  The  best  method  that  has  been  found  is  to  assume,  with 
Ogura,  that  dw/dx  *  0  at  x  *  0  (a  reasonable  assumption  for  axial  sym- 
metry)  whereby  it  is  easily  shown  that  w  *  d  ty^/dx  at  x  ■  0.  With 
this  formulation  the  magnitude  of  w  seems  to  increase  sharply,  by  com¬ 
parison  with  the  rectilinear  model,  over  the  last  grid  interval  as  the 
axis  is  approached,  but  not  so  sharply  as  with  other  axisymmetric  formu¬ 
lations.  The  maximum  updrafts  shown  in  Fig.  1  are  all  on  the  central 
axis,  but  the  nature  of  Fig.  1  would  not  be  changed  if  the  magnitude 
of  the  maximum  updraft  not  on  the  central  axis  had  been  plotted  instead. 
The  effect  at  26  minutes  would  be  to  reduce  the  value  from  11.9  m  sec 
to  7.6  m  sec  ^  in  the  axisymmetric  model  and  from  1.88  m  sec  ^  to  0.77  m 
sec"^-  in  the  rectilinear  model.  These  points  are  indicated  on  the  fig¬ 
ure  by  the  symbol  X. 

Ogura  did  not  discuss  the  actual  magnitude  of  the  vertical  currents 
but  only  the  ratio  of  maximum  downdraft  to  maximum  updraft.  These  ratios 
for  the  present  experiment  are  shown  in  Fig.  2.  As  Ogura  reported,  the 
ratio  for  the  rectilinear  model  is  substantially  larger  than  that  for 
the  axisymmetric  model,  but  the  ratios  are  not  at  all  constant.  Both 
show  a  marked  decrease  for  the  first  7  or  8  minutes,  followed  by  a 
marked  increase.  If  isotachs  of  the  vertical  component  of  velocity  or 
streamlines  are  plotted  for  each  time  step,  it  is  found  that  the  center 
in  which  the  maximum  downdraft  is  located  during  the  first  7  or  8  min¬ 
utes  is  not  the  same  as  that  for  later  times.  Initially  there  is  a 
single-cell  circulation  pattern  centered  near  x  *  1600  m,  with  updraft 
for  smaller  values  of  x  and  downdraft  for  larger  values.  As  the  simu¬ 
lated  cloud  develops,  this  cell  strengthens  and  its  center  moves  inward 
and  upward,  but  a  weak  secondary  center  develops  and  moves  outward.  At 
first  the  maximum  downdraft  is  outside  the  secondary  center,  but  gradu¬ 
ally  the  downdraft  between  the  centers  incr'^oes,  and  after  about  8  min¬ 
utes  it  becomes  the  dominant  downdraft.  Thereafter  the  maximum  updraft 
and  maximum  downdraft  are  closely  related  to  each  other  and  to  the  area 
of  condensation,  although  secondary  maxima  and  minima  occur  at  various 
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spots  throughout  the  computational  region.  Hence  after  about  10  min¬ 
utes  comparisons  should  be  significant. 

The  ratios  shown  in  Fig,  2  are  generally  comparable  to  the  ratios 
of  .25  and  .06  reported  by  Ogura  at  8  minutes,  but  they  tend  to  be  some¬ 
what  higher.  The  ratio  of  the  ratios  in  Ogura's  experiment  was  4:1, 
whereas  in  the  present  experiment  it  is  3:1  or  2:1.  These  differences, 
which  do  not  appear  to  be  significant,  can  be  ascribed  to  differences 
in  the  models,  or  more  particularly  to  differences  in  initial  data. 

What  is  important  is  that  the  ratio  of  maximum  downdraft  to  maximum  up¬ 
draft  is  larger  in  the  rectilinear  model  than  it  is  thought  to  be  for 
actual  clouds,  whereas  it  has  a  more  realistic  value  in  the  axlsymmetric 
model. 

Study  of  the  streamlines  and  isotachs  has  disclosed  no  apparent 
reason  for  the  regular  oscillation  of  the  upper  curve  of  Fig.  2  after 
8  minutes,  other  than  that  it  is  connected  in  some  way  with  the  gravity 
waves  of  period  near  10  minutes  that  are  typical  of  this  convective 
model.  The  waves  are  probably  the  consequence  of  the  assumption  of  a 
rigid  upper  boundary.  It  is  not  known  why  the  lower  curve  of  Fig.  2 
fails  to  show  the  same  oscillation. 

Plots  of  temperature  departure  as  a  function  of  position  at  vari- 

•ft 

ous  times  (e.g.,  Fig.  3  }  show  a  strong  maximum  in  the  region  of  con¬ 
densation,  surmounted  by  a  cap  of  temperature  deficit  resulting  from 
forced  lifting  of  yet -unsaturated  air  above  the  active  cell.  This  cap 
tends  to  extend  downward  in  a  narrow  band  along  the  edge  of  the  cloud, 
where  evaporation  is  taking  place.  Farther  out  there  is  another  maxi¬ 
mum,  strong  in  the  rectilinear  model  and  weak  in  the  axlsymmetric,  re¬ 
lated  to  subsidence  of  relatively  dry  air.  It  is  this  temperature  ex¬ 
cess  together  with  the  temperature  deficit  associated  with  dry  adiabatic 
ascent  below  the  cloud  that  cause  the  development  of  the  reverse-circu¬ 
lation  cell  mentioned  by  Ogura.  In  the  present  experiment  this  coll 

“ft 

Figures  3,  4,  5,  and  7  are  taken  from  computer-constructed  charts. 
On  them  small  squares  indicate  grid  points  of  maximum  value  relative  to 
adjacent  surrounding  grid  points,  and  small  circles  represent  relative 
minima  in  the  same  sense.  The  ordinate  is  height  and  the  abscissa  is 
radial  distance,  both  in  meters.  The  edge  of  each  square  block  is  five 
mesh  lengths. 
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appeared  both  with  the  axisymmetric  and  the  rectilinear  models,  but  it 
is  very  weak  in  both  (set  Fig.  u  ) .  It  anything,  it  is  a  little  stronger 
in  the  axisymmetr Ic  model,  but  so  are  all  other  circulations  *  A  strong 
reverse  cell  has  appeared  ir.  previous  runs  of  the  rectilinear  model, 
however,  when  the  initial  impulse  was  a  temperature  excess  as  in  Ogura's 
experiment.  It  is  believed  that  its  failure  to  appear  strongly  here  is 
related  to  the  type  of  impulse. 

The  center  of  excess  temperature  in  the  subsiding  air  has  another 
effect  that  has  been  evident  in  most  models  of  this  type.  As  the  model 
cloud  develops,  its  active  vortex  moves  upward,  leaving  a  relatively 
passive  mass  of  cloud  at  lower  levels.  The  inflow  associated  with  the 
vortex  brings  in  warm,  dry  air,  causing  that  portion  of  the  cloud  to 
evaporate.  There  frequently  results  a  rather  globular  cloud  ir.  the  re¬ 
gion  of  the  active  vortex,  having  a  narrow  stem  extending  down  along, 
the  central  axis.  Because  the  rectilinear  model  has  a  stronger  subsi¬ 
dence  temperature  excess,  this  effect  is  particularly  noticeable  in 
models  with  such  a  geometry;  see  Ogura's  Fig,  20.  The  same  action  oc¬ 
curs  in  the  axisymmetric  model,  but  it  is  less  intense.  The  effect  is 
illustrated  in  Fig.  5,  which  is  taken  from  a  case  with  J  more  extensive 
initial  impulse  than  that  of  Figs.  '  and  4,  and  consequently  with  more 
vigorous  development.  In  many  instances  tne  stem  of  the  cloud  in  the 
rectilinear  model  Is  o..iy  a  single  mesh  unit  wide. 

The  continued  maintenance  of  a  narrow  stem  below  the  active  pare 
of  the  cloud  brings  one  to  a  consideration  of  entrainment.  In  other 
types  of  cloud  models  an  explicit  entrainment  term  is  included,  which 
generally  has  the  effect  of  mixing  some  of  the  air  from  a  uniform  en¬ 
vironment  with  uniform  cloud  air,  resulting  in  a  cloud  with  new  uniform 
properties.  In  the  type  of  model  considered  here  the  cloud  is  not  an 
entity,  but  only  a  collection  of  grid  points  having  in  common  the  single 
property  r  >  0.  Hence,  the  ordinary  type  of  entrainment  term  is  not 
pertinent . 

There  are,  however,  two  types  of  entrainment  Ln  this  model.  Dy¬ 
namic  entrainment  refers  to  the  process  whereby  a  particle  trajectory 
leads  from  a  point  outside  the  cloud  to  a  point  within  it.  This  occurs 
when  a  parcel  without  liquid  water  reaches,  through  some  process,  a 
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state  of  saturation,  and  condensation  commences.  Such  a  parcel  has 
then  been  entrained.  Turbulent  entrainment  is  mere  comparable  to  the 
normal  notion  of  entrainment--a  form  of  mixing.  In  this  model  it  can 
occur  explicitly  through  the  Fickian  eddy  diffusion  terms  or  implicitly 
through  the  interpolation  required  by  the  Lagrangian  part  of  the  com¬ 
putation. 

Since  u  =  0  at  x  =  0,  dynamic  entrainment  cannot  occur  along  the 
central  axis  (except  trivially  at  the  bottom  and  top  of  the  cloud). 

It  can,  however,  occur  at  all  other  points  about  the  cloud,  and  it 
creates  the  mushroom  effect  noted  above.  Turbulent  mixing,  on  the 
other  hand,  can  occur  at  any  point  of  the  cloud,  including  the  central 
axis,  arid  it  may  be  the  key  to  control  of  unrestrained  growth  both  in 
the  model  and  in  nature.  The  basic  sounding  used  in  the  present  experi¬ 
ment,  the  lower  part  of  which  is  shown  in  Fig,  6,  is  only  slightly  less 
stable  and  more  humid  than  the  mean  August  Caribbean  sounding.  A  few 
cumulonimbus  clouds  were  observed  over  the  ocean  at  the  time  of  this 
sounding,  but  most  of  the  clouds  grew  only  to  moderate  heights  (Simp¬ 
son,  Simpson,  Andrews,  and  Eaton,  1965).  Yet  the  parcel  method,  which 
does  not  allow  for  entrainment,  would  predict  that  any  disturbance  would 
grow  to  the  tropopause .  It  is  apparent  that  nature  uses  entrainment 
to  check  runaway  convection. 

I';  in  the  present  model  the  eddy  exchange  coefficients  are  zero, 
the  development  along  the  central  axis  is  almost  identical  to  that  of 
the  parcel  method.  Computed  soundings  at  the  later  stages  of  develop¬ 
ment  show  an  almost  dry-adiabatic  lapse  rate  below  the  cloud  and  a  sat¬ 
urated-adiabatic  lapse  rate  within  the  cloud.  This,  of  course,  does 

not  hold  away  from  the  central  axis,  where  dynamic  entrainment  occurs. 

2  -1 

The  present  computer  runs  were  made  with  =  v^,  =  v  =  40  m  sec  , 

the  smallest  value  Ogura  found  to  have  any  discernible  effect  on  the 

■>v 

computations.  It  is  planned  to  repeat  the  computations  with  larger 


Molencarnp  (1968)  showed  that  with  Ogura's  finite-difference  scheme 
there  was  also  an  implicit  eddy  exchange  equivalent  to  the  explicit  ex¬ 
change  with  vM  “  40,  The  present  finite-difference  scheme  is  differ¬ 
ent,  but  the  effect  is  probably  comparable. 
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and  perhaps  variable  coefficients  to  see  whether  this  will  curb  the 
model's  tendency  toward  excessive  updrafts  and  temperature  departures 
within  the  cloud.  This  is  particularly  important  for  the  axisymmetric 
model,  for  its  strong  vert'^al  wind  along  the  axis  can  lead  to  compu¬ 
tational  instability  and  the  necessity  to  terminate  computation  before 
the  cloud  has  reached  its  mature  stage. 

Some  types  of  models,  such  as  that  of  Simpson,  Simpson,  Andrews, 
and  Eaton  (1965),  require  specification  of  cloud  diameter,  which  con¬ 
trols  the  rate  of  entrainment,  and  consequently  determines  the  height 
to  which  the  cloud  will  rise.  In  the  present  model,  by  contrast,  cloud 
diameter  is  one  of  the  results  that  come  out  of  the  computation.  It 
may  be  related  to  the  diameter  of  the  initial  impulse,  but  in  the  many 
runs  that  have  been  made  it  has  invariably  been  less  than  that  diameter. 
It  is  interesting  to  note  that  with  identical  initial  impulses  the  di¬ 
ameter  of  the  axisymmetric  cloud  turns  cut  to  be  a  little  greater  than 
that  of  the  rectilinear  cloud.  The  formulation  of  this  model,  however, 
does  not  ensure  that  there  will  be  lesser  entrainment  for  a  greater 
diameter,  except  in  that  the  wider  simulated  cloud  has  many  more  in¬ 
terior  grid  points,  but  only  a  few  more  boundary  grid  points.  Thus 
the  ratio  of  mass  entrained  to  total  mass  is  less,  but  this  is  not  a 
quantity  that  is  considered  by  the  model. 

Ogura's  observation  of  the  development  of  a  center  of  positive 
vorticity  in  connection  with  the  temperature  maximum  associated  with 
sc  sidence  seems  not  to  be  verified  in  the  present  experiment.  Instead, 
the  vorticity  in  both  models  app.-ars  to  be  organized  into  several  hori¬ 
zontal  layers  o’  like  sign  (in  the  lower  part  of  the  computational  re¬ 
gion  at  least),  the  vorticity  center  associated  with  the  active  cloud 
circulation  being  embedded  in  one  of  t’’em.  Tins  is  illustrated  in 
Fig.  7,  which,  incidentally,  has  Its  sign  convention  opposite  from 
Ogura's.  The  vorticity  chart  for  the  axisymmetric  model  shows  a  stronger 
principal  maximum  and  more  noise  than  that  for  the  other  model,  but  they 
are  essentially  similar. 


Ansynrrjetrtc  Rectilinear 
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V.  CONCLUSIONS 


Despite  certain  fundamental  differences  between  the  two  approaches 
most  notably  in  the  computation  of  temperature  and  moisture  content, 
Ogura's  model  and  *"he  present  c ;e  give  essentially  similar  results. 

In  particular,  Ogura’s  observations  on  the  differences  between  models 
with  ax i symmetric  and  rectilinear  symmetry  are  generally  borne  out  by 
the  present  experiment,  even  though  Ogura's  comparative  runs  dirfered 
from  each  other  in  more  than  the  basic  geometry.  Ogura,  however,  pre¬ 
sented  only  fragmentary  results  of  the  comparison,  whereas  the  present 
study  involves  a  somewhat  more  systematic  investigation  into  the  ef¬ 
fects  of  basic  geometry  on  the  model. 

The  principal  conclusion  reached  by  Ogura  is  that  the  downdraft 
is  relatively  stronger  a  the  rectilinear  model  than  in  the  axisymmetri 
This  is  true  also  of  the  present  study,  for  the  ratio  of  maximum  down- 
draft  to  maximum  updraft  is  about  .35  to  .40  for  the  rectilinear  model 
nd  near  .15  for  the  other.  These  ratios  are  not  constant  with  time, 
however,  and  their  physical  significance  seems  to  change  ab.  it  the 
eighth  minute.  At  later  times,  however,  the  comparisons  are  consistent 

What  is  even  more  striking  than  the  ditference  in  ratio  of  down- 
draft  to  updraft  is  the  difference  in  speed  and  strength  of  development 
In  the  present  experiment  the  ax isvmme t r ic  model  showed  much  stronger 
updrafts  than  the  rectilinear,  and  consequently  its  model  cioud  was 
both  bigger  and  wetter.  As  a  matter  of  fact,  its  growth  was  such  that 
the  updrafts  became  unrealistically  large,  requiring  termination  of 
computation  before  the  model  cloud  could  reach  its  mature  stage.  This 
was  partly  the  result  of  the  instability  and  high  humidity  of  die  ini¬ 
tial  sounding  used.  Experience  has  shown  that  a  well-placed  st.ble 
layer  can  easily  overcome  the  unrestrained  growth.  It  has  also  shown 
that  the  rate  oi  growth,  though  not  its  final  extent,  can  he  control. led 
by  varying  the  si. .2  of  tue  initial  impulse.  These  observations  apply 
mainly  to  the  axisyrsnetrlc  model;  the  much  slower  rectil ineu.  model 
was  not  run  long  enough  in  this  experiment  to  develop  unrealistic  up¬ 
drafts,  although  the  indications  wete  that  it  would  eventually  do  jest 
that,  Ogura  did  not  report  this  particular  difficulty,  but  his  compu¬ 
tational  region  was  so  small  Chat  the  model  cloud  impinged  on  the  upper 


boundary,  making  further  computation  unrewarding  from  the  standpoint 
of  realism.  The  implication  is  that  the  present  model  did  not  contain 
sufficient  turbulent  entrainment  to  slow  down  its  development.  It  ap¬ 
pears  that  through  turbulent  entrainment,  nature  is  able  to  prevent 
runaway  growth  of  convective  clouds  in  moist,  unstable  air.  This  ap¬ 
proach  will  be  investigated  in  the  course  of  future  work  with  the  model, 
but  the  strong  difference  between  the  two  geometries  remains,  and  is 
likely  to  continue  to  do  so. 

Aside  from  the  lack  of  dynamic  entrainment  on  the  central  axis, 
a  sort  of  feedback  mechanism  seems  to  ensure  the  faster  growth  of  the 
updraft  of  the  axisymmetric  model.  As  was  mentioned,  the  built-in  con¬ 
vergence  of  the  coordinate  system  ensures  that  the  updraft  will  initi¬ 
ally  be  relatively  stronger  than  the  downdraft.  When  condensation  oc¬ 
curs,  there  is  a  considerable  increase  of  temperature,  and  the  contri¬ 
bution  to  (9)  of  the  buoyancy  term  becomes  large  in  the  limited  region 
of  the  cloud.  This  leads  to  a  stronger  updraft,  more  condensation, 
and  so  on.  Thus  the  difference  due  merely  to  geometry  is  magnified. 
Between  the  cloud  edge  and  the  subsidence  temperature  maximum,  however, 
the  effect  of  the  buoyancy  term  is  negative,  but  it  is  not  enhanced  by 
condensation  or  evaporation.  Moreover,  the  convergence  of  the  coordi¬ 
nate  systems  is  less,  and  the  two  models  are  more  nearly  alike.  The^ 
result  is  that  even  though  the  rectilinear  model  shows  a  more  pronounced 
subsidence  temperature  maximum,  the  pattern  of  the  vortlclty  distribu¬ 
tion  is  not  much  different  in  the  two  models,  nor  is  the  pattern  of  the 
streamlines,  though  there  are  differences  in  the  spacing  of  the  stream¬ 
lines.  This  result,  however,  may  be  more  closely  connected  with  the 
initial  Impulse  than  with  the  geometry.  Ogura  noted  a  distinct  differ¬ 
ence  in  the  vorticity  patterns  and  the  development  of  a  reverse-flow 
cell  beneath  the  cloud  in  the  rectilinear  model.  The  present  model, 
when  run  with  a  temperature  impulse  similar  to  that  of  Ogura  rather 
than  one  of  humidity,  does  indeed  show  the  development  of  such  a  re¬ 
verse  cell  more  strongly  with  the  axisymmetric  geometry  than  with  the 
rectilinear.  On  the  whole,  though,  it  appears  that  the  humidity  im¬ 
pulse  leads  to  more  realistic  results,  so  it  has  been  retained. 
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Taken  together,  the  results  of  the  prese  i  experiment  coi  rm 
Ogura's  results,  but  go  beyond  them  to  suggest  that  the  two  geometries 
lead  to  strong  differences  in  the  simulated  clouds.  In  irost  respects 
the  axisynnetr ic  model  gave  results  more  in  keeping  with  obse-vations 
of  real  clouds,  but  it  still  needs  greater  restraint,  perhaps  through 
turbulent  entrainment,  to  prevent  runaway  growc'n.  Further  study  is 
suggested  in  that  respect  and  with  regard  to  the  effects  of  different 
initial  impulses  and  of  some  sort  of  precipitation  mechanism. 
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Appendix 


Let  cp  represent  any  of  the  primary  variables  (e.g.,  pressure). 
Then  by  definition 


cp  -  cp  +  '0B  +  cp ' 


where 

co  is  the  (constant)  mean  value  of  cp  over  the  entire  com¬ 
putational  region  at  initial  time  in  the  absence  of 
motion. 

cpB  “  cpg(z)  is  the  departure  of  cp  from  cp  at  a  given  height  z  at 
Initial  time  in  the  absence  of  motion, 
cp'  *  cp'(x,z,t)  is  the  departure  of  co  from  cp  +  cp_  at  any  location  (x,z) 

D 

and  any  time  t. 

The  symbols  used  in  this  report  have  the  following  meanings: 
Scalars 

B  ■  buoyancy  parameter 

L  ■  latent  heat  of  condensation 

R  ■  gas  constant  for  dry  air 

T  *  temperature 
* 

T  “  virtual  temperature 

^  ■  specific  heat  of  dry  air  at  constant  pressure 
Cpv  *  specific  heat  of  water  vapor  at  constant  pressure 
cw  *  specific  heat  of  liquid  water 
e&  ■  saturation  vapor  pressure 
g  ■  acceleration  due  to  gravity 
p  ■  pressure 

r^  ■  mixing  ratio  of  liquid  water  to  dry  air 
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r  =  satiation  mixing  ratio 
s 

=  mixing  ratio  of  water  vapor  to  dry  air 
t  =  time 

u  *  horizontal  component  of  wind 
w  -  vertical  component  of  wind 
x  radial  coordinate 
z  =  vertical  coordinate 
a  =  specific  volume 

e  =  ratio  of  molecular  weigh  -  of  water  vapor  and  dry  air 
T]  =  horizontal  component  of  vorticity 

=  coefficient  of  eddy  transport  of  momentum 
=  coefficient  of  eddy  diffusion  of  water  substance 
vT  -  coefficient  of  eddy  diffusion  of  temperature 
t  =  stream  function 

Vectors 

V  =  ~ind  velocity 
k  =  vertical  unit  vector 
ii'  =  vorticity 
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